Squeezing Out the Entropy of Fermions in Optical Lattices 
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At present, there is considerable interest in using atomic fermions in optical lattices to emulate 
the mathematical models that have been used to study strongly correlated electronic systems. Some 
of these models, such as the two dimensional fermion Hubbard model, are notoriously difficult to 
solve, and their key properties remain controversial despite decades of studies. It is hoped that 
the emulation experiments will shed light on some of these long standing problems. A successful 
emulation, however, requires reaching temperatures as low as 10~ 12 K and beyond, with entropy 
per particle far lower than what can be achieved today. Achieving such low entropy states is an 
essential step and a grand challenge of the whole emulation enterprise. In this paper, we point out 
a method to literally squeeze the entropy out from a Fermi gas into a surrounding Bose-Einstein 
condensed gas (BEC), which acts as a heat reservoir. This method allows one to reduce the entropy 
per particle of a lattice Fermi gas to a few percent of the lowest value obtainable today. 
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Currently, many laboratories are trying to realize the 
anti-ferromagnetic (AF) phase of the 3D Hubbard model 
using ultra-cold fermions in optical lattices Q. Although 
the AF phase is well known in condensed matter, its 
realization will be a major step for the emulation pro- 
gram since it requires overcoming the serious challenge 
mentioned above which is common to all cold atom 
emulations 0. To understand the origin of the prob- 
lem, recall that the strongly correlated states of lattice 
fermions emerge in the lowest Bloch band of the optical 
lattice. To put all the fermions in the lowest band, a 
sufficiently deep optical lattice is required. In such deep 
lattices, many known methods of cooling fail. For exam- 
ple, standard evaporative cooling does not work because 
the magnetic repulsive potential used in the evaporation 
process is not strong enough to overcome the deep lattice. 
As a result, all current experiments on lattice quantum 
gases resort to the conventional cooling scheme, first cool- 
ing the quantum gas in a harmonic trap (without optical 
lattice) to the lowest entropy state possible, and then 
turning on the lattice adiabatically[3][4]. The hope is 
that one could reach the state of interest when the lat- 
tice depth is sufficiently high. 

To realize these strongly correlated states with the 
current scheme, it is necessary that the entropy of the 
gas prior to switching on the lattice be less than that of 
the strongly correlated state one wishes to achieve. For 
3D fermion Hubbard models, recent studies [H show that 
it is possible to reach the antiferromagnetic (AF) phase 
slightly below the Neel temperature Tjv with the conven- 
tional cooling method. A similar calculation by Tremblay 
et.al.Q for the 2D Hubbard model, however, showed that 
the conventional scheme cannot reach even the pseudo- 
gap regime, which exists at a higher temperature than 
the anticipated superconducting phase. Although these 
calculations are for homogenous systems, they apply to 
confining traps as long as the majority of the sample is 
a Mott insulator with one fermion per site. These stud- 
ies show that even under optimal conditions, one can at 
most reach the AF phase close to the magnetic ordering, 



but not low enough temperatures to study ground state 
properties. 

The problem of the conventional method is that it can 
only cool atoms before they are loaded onto an optical 
lattice. There is no way to reduce the entropy further 
once the lattice is switched on. Thus, the lowest en- 
tropy attainable today in an optical lattice plus har- 
monic trap is just the lowest entropy achievable in a 
harmonic trap alone. In the case of Fermi gases, it is 
S/N = tt 2 (T/T f )[7}, where T F is the Fermi temperature 
in the trap. For T/Tp = 0.05, which is very much the 
limit today, we have S/N ~ 0.5. Any strongly corre- 
lated states with lower entropy are unreachable by this 
method. It is therefore important to find ways to reduce 
the entropy of the system significantly below this value. 



(A) OUR ENTROPY REMOVAL SCHEME: 

The purpose of this paper is to present a scheme to 
produce a lattice Fermi gas with about one fermion per 
site with entropy per particle much lower than what can 
be achieved today. Our method is based on the principle 
of entropy redistribution and the removal of entropy by 
isothermal compression. It consists of the following steps: 

(I) We immerse the lattice fermions in a BEC which acts 
as a heat reservoir. The fermions are confined in a har- 
monic trap and a strong optical lattice. The bosons sees 
a much weaker trapping lattice than the one that traps 
for the fermions. They are confined in a loose trap and 
cover the entire fermion system. The traps for bosons 
and fermions are species specific [8] so that they can be 
varied separately. (See below for discussions of these po- 
tentials.) 

(II) We compress the fermion harmonic trap adiabati- 
cally to turn the fermions at the center into a band insu- 
lator, which has two fermions per site and essentially zero 
entropy During this process, a substantial amount of the 
original fermion entropy is pushed into the bosons, while 
the entire system has little temperature change because 
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of the large heat capacity of the BEC compared to the 
lattice fermions. Hence, even though the process is adi- 
abatic for the entire Bose-Fermi system, it is essentially 
isothermal as far as the fermions are concerned. 

(III) After pushing out the fermion entropy into the 
BEC, we remove it by evaporating away the bosons all at 
once, leaving the remaining fermions to equilibrate. Since 
the band insulator is incompressible, only its density and 
entropy near the surface are affected during this process. 
This in turn severely limits entropy re-generation dur- 
ing the equilibration process. As a result, the entropy of 
the re-thermalized band insulator has a similar ultra-low 
value as before boson evaporation. 

(IV) We open up the the fermion harmonic trap adiabat- 
ically to lower the density of the lattice fermions. In this 
way, one can produce a Mott insulator or other states 
with fractional filling with the same ultra-low entropy. 

Before proceeding, we return to discuss the construc- 
tion of the trapping potentials mentioned in (I) . To have 
an optical lattice that confines the fermions tightly and 
the bosons loosely, the energy difference between the 
ground state (S'-state) and the excited state (P-state) 
of fermions must be smaller than that of the bosons, 
AEf < AEf,. In this way, one can choose a laser 
with frequency (oj) red detuned with respect to both 
excitation energies, (fuv < AEf, AEf,), such that the 
detuning for the fermions is smaller than that for the 
bosons, AEb — fiw > AEf — %u). For 40 K fermions, its 
45 to 4P transition has a wavelength 740nm, where as 
the wavelength of the 2S to 2P transition of 7 Li boson 
is 671nm, which satisfies the aforementioned condition. 
Moreover, the difference between these two excitation en- 
ergies are large enough so that the detuning AEf, — hu 
and AEf — Tilo sufficiently large to suppress heating due 
to spontaneous emission. Condition (I) can therefore be 
satisfied. 

Another scheme that makes use of the hyperfine struc- 
ture of the P state of 87 Rb was pointed out in ref.[||. 
By tuning the laser frequency to 790.01nm, which is be- 
tween the two hyperfine states P 3 / 2 and P\/i of 87 Rb, 
Rb bosons sees no lattice because the red detuned lattice 
due to P3/2 is cancelled by the blue detuned lattice due 
to Pi/2- On the other hand, such laser will generate an 
attractive potential for 6 Li and 40 K, since its wavelength 
is longer than those of the S-P transitions of 6 Li and 40 K 
respectively. 

To illustrate our scheme, we shall first discuss the basic 
properties of lattice fermions and the important process 
of entropy-redistribution. 



(B) NUMBER DENSITY AND ENTROPY 
DISTRIBUTIONS OF LATTICE FERMIONS: 

A Fermi gas in the lowest band of a (3D) optical lattice 
is described by the Hubbard model 

K = J + W - fj,N, (1) 

where J = — JJ2{n r/) a a R er a R' a describes hopping 
of fermions with spin a, (a =T,J,) between neighbour- 
ing sites R and R', J is the tunneling integral, W — 
^Er^ri^ri describes the on-site repulsion (U) be- 
tween spin up and spin down fermions; a R and n R a = 
a R rr a R a are the creation and number operators of a 
fermion with spin a at site R; N = J2n a n R a * s the 
total fermion number, and [i is the chemical potential. 

The possible states of the fermions are: band insulator 
(BI), Mott insulator (MI), and "conducting" state (C), 
corresponding to two, one and a non-integer number of 
fermions per site, respectively. The Mott insulator will 
develop AF order at the Neel temperature T/v ~ J 2 /U, 
while the conducting state is expected to have a su- 
pcrfluid ground state. The strongly correlated regime 
emerges when 

U » J» J 2 /U, (2) 

with J 2 /U being the smallest energy scale. At present, 
experiments operate in the temperature range 

U >T> J. (3) 

For simplicity, we set Boltzmann's constant ks = 1. Our 
goal is to perform operations in this high temperature 
regime so that the fermions will lose a substantial amount 
of entropy. 

Within the temperature range described by eq.Q, J 
can be treated as a perturbation and K is site-diagonal 
to zeroth order in J. The number occupation at site R 
is then 

2e' l/T + 2e { - 2 ^ u ^ T 
n R (T, n) = }Jn R , ff > = 1 + 2eM/T + e ( 2M -r/)/T - ( 4 ) 

(T 

The entropy per site is s R = d(T\nZ R )/dT, where 
Z n (T, fj,) = Tre-x*/ T = 1 + 2e^ T + e (^-u)/T is the 
partition function at site R. Explicitly, we have 

s r( T : M) = lnZ R (T, /x) + (E R (T, fi) - nn R (T, fi))/T (5) 

and P R (T, M ) = C/(n RT n Ri ) - e^- u V T /Z n (T, M ). 

In a harmonic trap V^(R) = Mcu 2 R 2 /2, the density 
and entropy distributions in the temperature range U > 
T > J can be calculated from eq.((4]) and eq.((5]) using 
local density approximation (LDA) by replacing [i with 
/i(R) = fj, — V^(R). A typical distribution is shown in 
Figure 1. The chemical potential fi and temperature T 
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FIG. 1: Density distribution np^ (red) and entropy distribu- 
tion sj^ (blue) of a Fermi gas in the lowest band of an optical 
lattice and in the presence of an harmonic trap. BI and MI 
denote band insulator and Mott insulator, respectively. The 
regions labeled CI and C2 are the "conducting" regions where 
the fermions have large number fluctuations and are mobile. 



T > J using eqs.(H|), ([5]), and ([6]), it is useful to under- 
stand them using simple estimates. The total entropy of 
the system is 

S = AnRlAR^ + 4irR 2 2 AR 2 J + ^{Rl - Rlf-f, 

(8) 

where d is the lattice spacing, and ARi and Ai? 2 are the 
widths of the (CI) and (C2) regions. They are given by 
A/i(i?i) - T and A[i(R 2 ) ~ T, or 

AR l =T/(Muj 2 R l ) AR 2 — T/(Mlu 2 R 2 ). (9) 

The quantities sc and sm are the average entropy per 
site in (C) and (MI) regions: ln3 < sc < ln4, and ln2 < 
~§m < ln4. From eq.®, we also have 

R\ - R\ = 2U/(Mu 2 ), (10) 

showing that the {MI) region shrinks as w increases. 

(C) ENTROPY LOCALIZATION AND 
REDUCTION: 



m cq. O and (|5| are determined from the number and 
entropy constraints, 

N = E "r( t - t* - K,(R)). S = s r( t > A* - V U (R)), 

R R 

(6) 

where N is the total number of fermions and S is the to- 
tal entropy produced in the convention cooling scheme 
before the lattice is turned on. Figure 1 shows the 
following phases: [(BI) : n R — ► 2, for U < /i(R)]; 
[(MI) : n R -> 1 for < fi(R) < U}; [(Vacuum) : n R -> 
for fi(R) < 0]. There are also two "conducting" phases 
(CI) and (C2) with non-integer number of fermions per 
site. They are (CI) : 1 < n R < 2, with /x(R) ~ [/, and 
(C2) : < n R < 1 with fi(K) - 0. We shall denote the 
centres of the conducting regions (CI) and (C2) as R l 
and R 2 . They are determined by 

H - \^Mu?R\ = U and fi - ^Muj 2 R 2 2 = 0. (7) 

The (BI) region has essentially zero entropy per par- 
ticle, sr ~ 0, since the pair (||) is the only possi- 
ble configuration at each lattice site. At temperatures 
U > T > T N , (MI) has spin entropy sr ~ ln2, since 
both f and |, are equally probable. At the same temper- 
ature range, (CI) has sr ~ ln3, since the doublet 
and the single spin states, f and j are all equally proba- 
ble. At even higher temperatures, sr of both (MI) and 
(CI) can rise as high as ln4, as the probability of hav- 
ing an empty site increases from zero. Similar situation 
occurs in (C2). 

While the total entropy and other properties of the 
system can be calculated in the temperature range U > 



From Figure 1, we see that if we compress the trap 
adiabatically (by increasing lo), the band insulator will 
grow at the expense of other phases. As a result, the 
Mott regions and the "conducting" regions, and hence 
all the entropy, are pushed to the surface, as shown in 
Figures 2a and 2b. Since all the entropy in the bulk is 
squeezed into the surface layer, the entropy density at 
the surface must rise above its former value to keep the 
total entropy constant, which means the temperature the 
system will increase, as shown in Figure 2. The fact that 
adiabatic compression causes heating is a consequence of 
thermodynamics [2| . 

As compression proceeds, it reaches the point where 
the width of the (MI) region becomes so thin that it is 
of the order of a few lattice spacings. (See Fig. 2a and 
2b). This occurs at R 2 — R\ « R\ ~ R 2 ~ R, where 
R is the radius of the (BI) region related to the total 
fermion number N by N = 2(47r/3) (R/d) 3 . In this case, 
we have 

N 2\ZN) \ Mu 2 d 2 /' 1 ' 

Eq.|jni) shows that T rises as w during an adiabatic 
compression. 

If, however, we are able to keep the temperature con- 
stant during the compression, the entropy density at the 
surface will remain at its initial value, (see Figure 2c, 2d). 
As a result, the total entropy will be reduced by the same 
amount as the reduction of entropic volume (i.e. the re- 
gion that contains entropy), which is substantial as the 
latter changes from a 3D volume to a surface layer. This 
can also be seen in eq. (fTTj) . where S/N drops as lq~ 2 in 
an isothermal compression. 
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As isothermal compression continues, the widths of the 
(C) and (MI) regions will shrink down to a lattice spac- 
ing d, and the entropic surface region will eventually be 
reduced to a surface of thickness d. At this point, LDA 
breaks down. Due to thermal fluctuation, the thickness 
of this surface layer will remain approximately one lat- 
tice spacing, regardless of further increase of the trapping 
potential. As long as T is above the spin ordering tem- 
perature T N , the entropy per site is still ln2. This limits 
the lowest entropy attainable in isothermal compression 
for T > J 2 /U to S* = (ATrR 2 d)(s/d 3 ), or 

S* 3 / 8tt \ 1/3 

— = - s, (12) 

where s ~ ln2. For N = 5 x 10 6 (3 x 10 7 ), S*/N is 0.012 
(0.007) which is 2% (1%) of the lowest value attainable 
today. 

(D) ENTROPY TRANSFER BETWEEN 
FERMIONS AND BEC: 

To keep the fermions at roughly constant temperature, 
they must be in contact with a heat reservoir. A BEC 
is an ideal medium for this purpose, for it has a much 
higher heat capacity than the fermions and can therefore 
keep the temperature roughly constant. In this section, 
we shall provide explicit calculations to demonstrate this 
fact. Our results show that the transfer of entropy be- 
tween fermions and bosons is accurately described by the 
simple isothermal compression model in Section (C). 

We shall consider a mixture of BEC and fermions. 
The BEC is contained in a harmonic trap V B (R) = 
^M B u B Il 2 with frequency uj b , where M B is the boson 
mass, and another trap V F (R) = \M f lo 2 f T{? is used for 
the fermions. The fermions are also confined in an optical 
lattice. We shall assume these potentials are species spe- 
cific as described in summary (I) in Section (A) , so that 
V F can be varied independently of V B . This is important 
because the compression of fermions should not lead to a 
substantial compression of the BEC. Otherwise the tem- 
perature of the BEC will rise, making it less efficient in 
absorbing the entropy of the fermions. 

To study the entropy transfer between fermions and 
bosons, we first consider a homogenous Bose-Fermi sys- 
tem. In the grand canonical ensemble, the hamiltonian 
is 

K = Kf(hf) + K b {hb)+H bf , (13) 

where Kp is the Hubbard hamiltonian eq. JTJ) with [i now 
denoted as [ip. K B is the hamiltonian for bosons with 
chemical potential /i B , 

k b = j (^W 1 -vi + ^P&U-iiBpi), 

(14) 
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FIG. 2: Figure (2a) shows the number density nj^ of N = 
2 x 10 6 40 K fermions in an optical lattice with lattice height 
V = 15Er and U = 0.5Er at harmonic trap frequencies 
v F = uj f /2tt = 100Hz (blue), 150Hz (purple) and 600Hz 
(red), calculated from eqs.Q-©. The entropy per particle 
is fixed at S/N = 0.5. The corresponding entropy distribu- 
tions sj> are shown in Figure (2b), with temperatures 27nK 
(blue), 88nK (purple), and 1880nK (red), respectively. As the 
trap frequency is increased from 100Hz (blue) to 150Hz (pur- 
ple), Figs. (2a) and (2b) show that the Mott phase melts away 
and the entropy density at the center of the conducting layer 
increases, leading to a rise in temperature. This rise is very 
rapid. It scales as uj 2 as seen in eq. (|lip . Figures (2c) and (2d) 
show the number density nj^ and entropy density of sp^ for 
the same parameters as Figure (2a) and (2b), but now with 
the temperature fixed (T = 50nK) instead of entropy per par- 
ticle. As the harmonic trap is compressed from 100Hz (blue) 
to 150Hz (purple) and then 600Hz (red), S/N decreases from 
0.732 (blue) to 0.333 (purple) and finally, 0.021 (red), reach- 
ing a value which is only 4% of the lowest value attainable 
today, and has not yet reached the limit shown in eq. (|12[) . 

where <jy is the creation operator for bosons and gsB = 
Airh 2 a B B / M B . Hbf is the boson-fermion interaction, 

a R 

with coupling constant g BF , d is the lattice spacing, and 
n FR (n B R ) is the number of fermions (bosons) in a unit 
cell centered at site R. 

To calculate the properties the system, we make the 
following approximations : 

(i) a mean field decomposition of H B p, replacing it by 
h bf = 9bf J2R(^F,R n B+ n Ff l B^ n F n B)^ where n F = 

V*F,R/) n B — \^B,R/) 

(ii) treating the hopping term J as a perturbation of K F 
as in Section (B). This is justified in the temperature 
regime U > T > J; 

(iii) applying the Hartree-Fock approximation for ther- 
modynamics of the bosonsQ. With approximation (i), 
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K becomes 

K' = k F (fj, F -g BF n B ) + k B (fj, B -g BF n F )-g BF n B n F . 

(16) 

The pressure P(T, fj,p, fx B ) = Q^ThaTre^ l T is then 
P(T,/j, B ,fi F ) = P F (T,n F - g BF n B ) 

+ P b( T ^B - 9BF n F) + 9BF n F n B,( 17 ) 

where P FB (T,^ FB ) = fT 1 Tlne~ i ^ ( ^ )/T , and Q is 
the volume of the system. The fermion and boson densi- 
ties are given by 



dP dP 
nF = , n B 



dfii 



(18) 



Within approximation (ii), the fermion density is readily 
given by eq. (J4J) as 



n F = n F {T,[i F - g BF n B ). 



(19) 



To apply approximation (iii), we follow the procedure in 
reference [9j to calculate the boson density for both the 
normal and superfluid parts of the boson cloud, which is 
of the form 



n B = n B (T,fj, B - g BF n F ) 



(20) 



according to eq. (fTrl)) . (The presence of the superfluid will 
depend on the value [i B — g BF n F ). 

Eqs. (fT9"|) to (|2"D|) form a complete set of equations that 
determine (n Bl n F ) self consistently as a function of 
(T, fi Bl /z F ). All these have to be evaluated numerically. 
These solutions then allow one to evaluate the entropy 
density which is s(T, fi F , [i B ) = or 



dP 
dT 



>BFi 



(21) 



where s F , s B , s BF are defined as fermion, boson, and "in- 
teraction" entropy density respectively, 



dP F 
dT 



dP B 
dT 



^BF 



9BF d -^, (22) 



which will also be functions (T, \x Bl fi F ). To find the den- 
sity and entropy distribution in the presence of the boson 
and fermion traps, we can use LDA to replace [i FB by 

In our numerical calculations, we consider 5 x 10 6 of 
40 K fermions in a lattice immersed in a BEC of 87 Rb with 
N B = 5 x 10 7 bosons. The frequency of the fermion trap 
is initially set at v F = uj/(2tt) = 100Hz. The bosons are 
confined in a loose trap with frequency v B — lo b / (2tt) = 
10Hz, so that it covers the entire fermion system. We 
take a BB = 5.45 nm, and g BF = g BB /2. The lattice 
height is V — 15En and the fermions have a Hubbard 
interaction U = 1.5Er, where En is the recoil energy. To 
demonstrate the transfer of entropy from the fermions to 
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FIG. 3: Density distribution n R and entropy distribution s R 
of bosons (blue) and fermions (red) in a mixture of BEC and 
lattice Fermi gas. The total number of fermions N F (Bosons 
Nb) is 5 x 10 6 (5 x 10 7 ). The trapping frequency of bosons 
is fixed at vb = 10Hz. The number density of bosons and 
fermions at fermion harmonic trap frequencies 100Hz, 200Hz, 
and 600Hz are shown in Figs.(3a)-(3c). The corresponding 
entropy densities are shown in Figs. (3d) to (3f). As the fre- 
quencies v F of the fermion harmonic trap increases, the edge 
of the band insulator sharpens and the entropy distribution 
becomes more and more narrowed, yet the temperature of 
the system changes very little as the fermion trap is tight- 
ened, from 60nK to 64.2nif, and to 65.6nK. As a result, the 
height of the entropy distribution remains roughly constant. 
The narrowing of the entropy distribution, however, leads to 
a rapid drop in the entropy per particle, from 0.691 to 0.184, 
to 0.016, as shown in Figs.(3d)-(3f). 



the BEC, we consider an initial state with initial tem- 
perature 60nK, corresponding to the entropy per particle 
S F /N F = 0.691 for the fermions. For these parameters, 
most of the bulk is already a band insulator, similar to 
that in Figure (2a) and the entropy density is accumu- 
lated at the surface. (See Fig. (3a)). 

When the fermion trap frequency v F increases from 
100Hz to 600Hz, we see from Figs. (3a)- (3c) that the edge 
of the band insulator becomes sharper, while the entropy 
density becomes concentrated at the surface (Fig. (3d)- 
(3f)). Our calculation also show that the fermion entropy 
S F = f Sj decreases whereas the boson entropy S B — 
J s B increases, while Sb F = J s BF remains much smaller 
than S F and Sb during the compression, and the entropy 
of the entire Bose- Fermi mixture is a constant. 

Due to the large heat capacity of the bosons, the over- 
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all temperature T rises only moderately, from 60nK to 
65.6 nK. This shows that the simple isothermal compres- 
sion model for the lattice fermions in Section (C) is a 
reasonable approximation of the entropy transfer process 
between the bosons and the fermions. At ^ F =600Hz, the 
fermion entropy per particle is S F /N F ~ 0.016, which is 
about 3% of the best estimate achievable with conven- 
tional methods. 



by N F = J2u n F( T{f) ~ K,(R))- These two rela- 
tions uniquely determine the final temperature and 
chemical potential fx F of the band insulator, from which 
one can calculate the final entropy. With the parameters 
we mentioned before, we find that the entropy increase 
in this process is negligible, only about a few percent of 
the value before evaporation; hence our conclusions sum- 
marized in Section (A). 



(E): EQUILIBRATION OF THE LATTICE 
FERMIONS AFTER BOSON EVAPORATION: 

After transferring the entropy of the fermions to the 
BEC as shown in Fig. (3c) at the temperature 65.6 nK(= 
T") and at trap frequency v F = 600Hz, we evaporate 
all the bosons suddenly so as to obtain a pure fermion 
system. The lattice Fermi gas will then relax to its equi- 
librium state corresponding to the new trap frequency 
lu f = 2ttv f . In this process, entropy will be generated 
and temperature will rise. However, because the band 
insulator is incompressible, and because the trap has a 
larger value of lo f , entropy generation is limited. We 
find (in the calculation below) that the increases in en- 
tropy and temperature are about a few percent of their 
values prior to evaporation. In other words, the entropy 
of the band insulator after it reaches equilibrium is es- 
sentially the same as that before boson evaporation, i.e. 
S F /Np ~ 0.02. This will be the entropy inherited by the 
Mott insulator that emerges from the band insulator as 
one decompresses the trap. 

The calculation for the above processes is as follows. 
Immediately after evaporation, the lattice Fermi gas has 
total energy 

E {1) = £ R e(T« M F -14(R)-g BF n s (R)) 

+ ErK,(R)MR), (23) 

where e(T, fj,) = U{hn^hn : i) is the internal energy per 
site of a homogeneous lattice fermion with hamiltonian 
K F at the temperature range U > T > J, and 

Ue (2»-U)/T 
< T > A = 1 + e ,/T + e (2^U)/T - ( 24 ) 

n B F (R) are boson and fermion densities prior to bo- 
son evaporation, which were calculated from eqs. (fT9|) and 
(f20|) . When the system finally relaxes to equilibrium, it 

will have a different temperature T^' and chemical po- 

( f) 

tential £i F . The energy of the final system is 

E U) = Er^,^-^(R)) 

+ Er V u (R)n F (T(/) ,rfP-V u (R). (25) 

Since energy is conserved in this process, we have EW = 
EW . In addition, the total number of particles is given 



FINAL REMARKS: 

We have introduced a method that allows extraction of 
a substantial fraction of the entropy of a Fermi gas in an 
optical lattice after a strong lattice is switched on. More- 
over, the extraction process is conducted at the temper- 
ature regime T > J, much higher than the Neel temper- 
ature TV ~ J 2 /U. While our method makes explicit use 
of the band insulator, it is applicable to any system which 
has an equilibrium phase with a large gap. The idea is to 
use the gapful phase to push away all the entropy in the 
bulk into a surrounding Bosc-Einstein condensed gas. It 
should also be stressed that although the final stage of 
our method involves evaporating away the Bose-Einstein 
condensed gas, it is different from the usual sympathetic 
cooling both in purpose and in function. First, our evap- 
oration must be preceded by the localization of fermion 
entropy in order to be effective. Second, the purpose 
of removing the bosons is not to decrease the energy of 
the system (as in the usual sympathetic cooling), but 
to remove the entropy it absorbed from the fermions. A 
natural question is whether one can efficiently reduce the 
entropy of a Mott phase which is in direct contact with 
a BEC by evaporating on the latter. The answer is neg- 
ative, as shown by our calculations. The reason is that 
the interactions between bosons and fermions typically 
have very weak spin dependence. The removal of bosons 
therefore has little effect on the spin entropy of the Mott 
phase, which will remain close to its value before boson 
evaporation, ln2 per particle. 

Finally, we point out that our scheme assumes that 
it is possible to expand the trap adiabatically to turn a 
band insulator into a Mott insulator. This implies that 
the relaxation time for particle re-distribution is suffi- 
ciently fast. Since mass transport occurs at the interface 
between different phases, (i.e. regions (CI) and (C2)), it 
will take place within the time scale Ti/J. It is therefore 
helpful to use lattices with sufficiently large tunneling, 
(say, around 10 — 15Er), and to reach the large U limit by 
increasing the interaction using a Feshbach resonance [10(. 
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